Linear and quadratic temperature dependence 
of electronic specific heat for cuprates 
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We model cuprate superconductors as an infinite layered lattice structure which contains a fluid 
of paired and unpaired fermions. Paired fermions, which are the superconducting carriers, are 
considered as noninteracting zero spin bosons with a linear energy-momentum dispersion relation, 
which coexist with the unpaired fermions in a series of almost two dimensional slabs stacked in their 
perpendicular direction. The inter-slab penetrable planes are simulated by a Dirac comb potential 
in the direction in which the slabs are stacked, while paired and unpaired electrons (or holes) are free 
to move parallel to the planes. Paired fermions condense at a BEC critical temperature at which 
a jump in their specific heat is exhibited, whose values are assumed equal to the superconducting 
critical temperature and the specific heat jump experimentally reported for YBaCu07_a; to fix our 
model parameters: the plane impenetrability and the fraction of superconducting charge carrier. 
We straightforwardly obtain, near and under the superconducting temperature T c , the linear (7 e T) 
and the quadratic (aT 2 ) electronic specific heat terms, with 7 e and a of the order of the latest 
experimental values reported. After calculating the lattice specific heat (phonons) C; from the 
phonon spectrum data obtained from inelastic neutron scattering experiments, and added to the 
electronic (paired plus unpaired) C e component, we qualitatively reproduce the total specific heat 
below T c , whose curve lies close to the experimental one, reproducing its exact value at T c . 



PACS numbers: 74.20. De, 74.25. Bt, 74.72.-h 



I. INTRODUCTION 

From the discovery of cuprate High Temperature Su- 
perconductors in 1986 (HTSC) many efforts have been 
made to explain the nature of their microscopic behavior 
as they are not completely described by the BCS theory 
0. The HTSC cuprates are perhaps the most studied 
both experimentally and theoretically until FeAs showed 
up. Their main characteristics can be summarized as 
follows: a small coherence length, usually of one or two 
nanometers; they have preferably either tetragonal or or- 
thorombic crystallographic structures; and the Cooper 
pairs, responsible for the superconductivity, move in the 
copper oxide planes which resemble a quasi-2D layered 
system. They also modify their structure Q and their 
critical temperature T c by changing the oxygen concen- 
tration, being this feature responsible for achieving or 
not superconductivity for a same compound with differ- 
ent oxygen concentrations. 

Particularly, the specific heat of the YBa2Cu307_ 2; 
cuprates, where x represents the oxygen dopage made 
with holes, has been widely studied and we want to point 
out four key characteristics which have been observed. 
First of all, even though it is not easy to observe, below T c 
there is a linear term j e T in the electronic specific heat, 
with 7 e the electronic specific heat parameter or Som- 
merfeld constant - sometimes referred to as 7(0)-, whose 
latest reported values are between 2 — 3 mJ/mol K 2 
This term is currently believed to come from the normal 
state electronic specific heat (see Ref. @ and references 
therein). In the second place, also at temperatures below 



T c , there is an aT 2 term in zero magnetic field 0, @, 0] 
(initially denied in some reports @), which changes to 
a H X / 2 T component in the presence of an external mag- 
netic field H, and is attributable to the superconducting 
part of the electronic specific heat. The reported values 
for this constant are of the order of tenths of mJ/mol 
K 3 [J], but as it happens for the linear term coefficient, 
the obtained values depend strongly on the conditions of 
each experiment and on the theoretical method each au- 
thor uses to subtract the other specific heat components. 
In the third place we have a "jump" (at zero magnetic 
field) in the specific heat at T c Q indicating a second 
order phase transition (which becomes a "peak" at finite 
magnetic field) and is widely believed to be the result 
of the influence of the superconducting electronic spe- 
cific heat C es . Finally, the total specific heat divided by 
temperature C/T shows an "upturn" (or "fishtail") for 
temperatures under 5 K, which several authors [5j claim 
has an intrinsic magnetic origin, present even at zero ex- 
ternal magnetic field. This latter low temperature fea- 
ture has been modeled as a contribution of five terms: 
normal electronic C en , superconducting electronic C es , 
lattice Ci , magnetic and hyperfine specific heats [f| [lC| , 
even though the dynamic mechanism beneath them is 
not completely known. 

Contrary to what one might think, to our knowledge, 
the lattice specific heat C\ of a cuprate is not fully de- 
scribed by standard theories. This specific heat is gener- 
ally considered as a contribution that doesn't change with 
the superconductivity onset and it has been observed to 
behave as a T 3 term below 5 K For the lattice spe- 
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cine heat a number of proposals have been made, such 
as: use of the Born- von Karman formalism (Tl| applied to 
the multiatomic anisotropic lattice; fits of Debye and/or 
Einstein models [l2l - fl4| ; the use of different polynomials 
or power series on a number of different variables 

EMI; 

adaptation of models that have given successful results 
in other elements (such as 4 He) |19l— : the use of the 
lower Landau level (LLL) formalism [2l| . to mention a 
few. 

On the other hand, experimentalists have been using 
indirect methods for obtaining the electronic specific heat 
C e - although a distinction between the normal and the 
superconducting part is sometimes not recognized - usu- 
ally by separating the lattice specific heat from the total. 
For example, Loram et al. [22j, construct the lattice spe- 



cific heat using a non-superconducting reference sample, 
using either a small variation on the oxygen content or in- 
troducing another element; or Meingast et al. (23|, who 
construct their phonon density of states usin g a local- 
density approximation; or Bessergenev et al. [l8j |. who 
develop a power series in terms of characteristic temper- 
atures which depend on phononic moments; or using an 
estimated phonon spectrum based on known lattice vi- 
bration frequencies and inserting it in a set of Einstein 
functions with characteristic temperatures (24|. In all 
these studies, the authors subtract the "lattice" specific 
heat obtained from the total experimental specific heat of 
the cuprate, and what is left is reported as the "electronic 
specific heat" , which appears as a small contribution, re- 
stricted mainly to the height of the jump of only a few 
percent (1 — 2%), and unwillingly transferring the intrin- 
sic uncertainties of their method to C e . Based on our 
results, we propose that this manipulation should be re- 
considered, since in this work we show that the electronic 
specific heat (normal and superconducting) contributes 
with a (30 - 40)% of the total. 

For conventional superconductors 7 e can be obtained 
from the linear term of the electronic specific heat in the 
normal state C en , while the superconducting component 
C es has an exponential behavior at very low temperatures 
predicted by BCS. For HTSC this term has been con- 
firmed to exist even in the superconducting state, how- 
ever, it is very difficult to separate it from the other com- 
ponents due to thermal fluctuations and to the need of 
very high external magnetic fields to suppress the charac- 
teristic upturn H| . The importance of determining the 
value of 7 e lies on its direct relation to the electronic den- 
sity of states and on the belief that it gives information 
about the interaction electron-phonon [5J. In addition, 
there are models, such as Anderson's Resonating Valence 
Bond (RVB) (2f|, that predict a linear term in the spe- 
cific heat, but until now, the controversy over whether 
7e comes from a residual zero-field term, from the su- 
perconducting electronic specific heat [1| or from normal 
electrons that didn't participate in the superconducting 
state (see, for example, [26| p. 85) goes on. In this work 
we adopt this point of view and consider the linear part 
of the electronic specific heat, j e T, is a result of unpaired 



electrons inside a layered structure, as we will show. 

Even though conventional superconductors do not ex- 
hibit an aT 2 electronic specific heat term, it has experi- 
mentally been confirmed that HTSC do [3| . Most of the 
reported values for this constant were obtained by first 
fixing 7 e and then fitting curves with other parameters 
[(| 0] , while others [J, [27j set an arrange of different exter- 
nal magnetic fields to cancel the interfering components. 
In the model presented in this paper, such a quadratic 
in temperature term in the electronic specific heat comes 
from paired electrons (superconducting state) , and values 
of a are of the order of some ones reported experimentally 

At zero external magnetic field, the "jump" AC in the 
specific heat at the transition temperature has also been 
reported with a great variety of values depending on each 
experiment. Currently, this feature is attributed to the 
electronic specific heat, its magnitude is of AC « 5 J/mol 
K according to several authors [3, EH Gjl H3| , and it has 
also been shown that for the same sample its jump di- 
minishes as the magnitude of an external magnetic field 
is increased @. As experiments have become more ac- 
curate, the shape of the jump has become sharper j4j. 
However, a roundness of the curves at the jump with 
positive curvature is usually justified by the finiteness of 
the samples and by the presence of thermal fluctuations, 
which are important in HTSC 28 1. 

In the framework of the most basic Boson-Fcrmion 
model [29M32| | of superconductivity, we assume Cooper 
pairs are composite-spin-zero-bosons with either zero or 
nonzero moments of center of mass, coexisting with a 
fermion fluid formed by the unpaired electrons. These 
Cooper pairs are pre-formed at some temperature T* > 
T c and can undergo a Bose-Einstein condensation (BEC) 
as temperature is lowered [3l], HH . We are aware that the 
number of preformed pairs increases as the temperature 
is lowered until T c , where its density is large enough to 
achieve coherence, independently of the mechanism by 
which the pairs are formed. Below T c we assume that 
the number of pairs remains constant. 

On the other hand, in order to include the effect of 
the layered structure of cuprates in the Boson-Fcrmion 
model, we previously calculated the BEC critical temper- 
ature and the thermodynamic properties for a system of 
non-interacting bosons immersed in a periodic multilayer 
array which represents the confinement agent [34], |35| . 
The multilayer array is simulated by an external Kronig- 
Penney (KP) potential at the delta limit case along the 
perpendicular direction to the Cu02 planes, while the 
particles are allowed to move freely in the parallel di- 
rections with an energy-momentum dispersion relation 
where the linear term predominates, as has been shown 
in [Hj]. 

Our system model begins with N electrons of mass 
m e interacting via a BCS type potential. There is a sub- 
group of electrons able to form pairs (Cooper pairs), since 
they are within a shell of width 2Tiujo around the Fermi 
energy Ep, where Tiujd is the Debye energy, coexisting 
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with a non-pairable group of electrons formed by those 
under and above the pairing shell, and are not eligible for 
pairing. From the first group, which we call the pairable 
electrons, we consider that only a fraction of them are 
paired, these are equal to a smaller fraction fN/2 which 
participate in the superconductivity; our assumption is 
based on the analysis of Uemura's plot (Fig. 2 of Ref 
[36| ) that shows that critical temperatures for cuprates 
are in the empirical range of T c w (0.01 - 0.06)T F [13 ■ 
With all the stated above in mind, the N electrons are 
grouped in three major components: paired electrons 
(boson gas) formed by a fraction / of half the total N 
electrons (inside the pairing shell) ; a fermion gas formed 
by the pairable but unpaired electrons (also inside the 
pairing shell); and the unpairable electrons (outside the 
pairing shell); plus a phonon gas due to the lattice. In 
Sec. Ill Al we obtain the grand potential from where it is 
possible to derive all the thermodynamic properties. This 
model depends on three physical properties: the sepa- 
ration between planes a; the impenetrability Pq of the 
planes, which is responsible for the anisotropy observed 
in cuprates; and the density of superconducting carriers 
/n/2, with n the fermionic number density. In section 
III B I we fix a with the experimental values reported and 
calculate the critical temperature of the boson gas made 
of Cooper pairs as a function of Pq and /. 



In section IIII Al we obtain the superconducting elec- 
tronic specific heat for the Cooper pairs fixing the un- 
known parameters Pq and / with the experimental T cexp 
and the magnitude of the "jump" in the electronic spe- 
cific heat at T c . As a consequence C es shows the observed 
T 2 behavior. Meanwhile, in section IIII Bl we derive the 
expressions for the normal electronic specific heat of the 
total unpaired (1 — f)N electrons (fermions) in a peri- 
odic layered structure, which shows the expected linear 
dependence on T. In Sec. IIII CI we add the two contri- 
butions and show that the total electronic specific heat 
at and under T c has the same behavior obtained by ex- 
periments. The electronic specific constants we obtain 
are 7 e = 5.2 mJ/mol K 2 and a = 4.3 mJ/mol K 3 , com- 
pared to 2.19 and 0.21 reported in Q, and to 25.1 and 
3.4 reported in [ijj]. In addition, since the lattice spe- 
cific heat represents the main contribution of the total 
specific heat, in section llVl we calculate the specific heat 
for the phonons using two different approaches: Debye 
model and phonon density of states from inelastic neu- 
tron scattering (INS) experiments (3^]. In section Ivl we 
sum these three specific heat contributions and compare 
the result with experiments. We find an excellent qual- 
itative agreement with the experimental C, where the 
contribution of the electronic specific heat is a significant 
part of the total. Finally, in Sec. IVII we present our 
conclusions. 



II. COOPER PAIRS IN A LAYERED 
STRUCTURE 

We begin by taking a group of pairable electrons im- 
mersed in a periodic layered array along the z direc- 
tion and free to move in the other two directions with 
a linear dispersion relation. The wave vector for the 
center of mass of the pair (CMM) is given by K = 
(K x ,K y ,K z ) ee ki + k 2 , while k ee l/2(k x - k 2 ) is the 
relative momentum, where ki and k 2 are the wave vec- 
tors of each electron of the pair. The solution for the 
Schrodinger equation for the pairs may be separated in 
the x — y and z directions, so the energy for each particle 
is Sk = £x x y + £x z ■ Here the total pair energy in the 
plane is sk x y = 2Ep — Ak , with Ak the binding energy 
of the pair for any temperature and any center of mass 
momentum. It has been shown that when K is non zero, 
but small, one can expand the binding energy from the 
Cooper equation in a series of powers [12], so the total 
energy in the plane is 

ex.,, = e + C x (Kl + if 2 ) 1 / 2 + 0(K 2 ), (1) 

where eo = 2Ep — Aq is a constant, C\ — {2 / it)h\i F2D is 
the linear term coefficient in 2D, vf2D is the correspond- 
ing Fermi velocity, Ao = 2hujrj exp(— 1/A) is the energy 
gap for K = and weak coupling (corresponding to the 
BCS theory), and A ee g(Ep)V the dimensionless cou- 
pling constant in terms of the electronic density of states 
at the Fermi sea g{Ep) and the non-local interaction V. 

Along the z-direction we use the Kronig-Penney [3^| 
potential, where the energies are implicitly obtained, as 
a function of a, from the transcendental equation 

Po(a/Ao) sin(aK z a)/aK z a + cos(a^a) = cos(fC z a), (2) 

with a 2 K ee 2me^/?i 2 , m = 2m e the mass of the 
composite-boson, and we have defined the dimensionless 
plane impenetrability Pq ee PX /a = mAXo/h 2 . The con- 
stant Ao = hj \/2irmkBTo is the de Broglie thermal wave- 
length of an ideal boson gas in an infinite box at the BEC 
critical temperature To = 27rK 2 n 2 3 / ' 3 /m/cB£(3/2) 2 / 3 ~ 

3.31h 2 n 2 J 3 /mks, with ns = N/(2L 3 ) the boson number 
density and A is the strength of the KP delta potentials 

Note that when Pq — > the energy goes to the free- 
particle energy sk z — > h 2 K 2 /2m in the z direction. Also, 
when we have small energies, ek z < h 2 iT 2 /2ma 2 , we get 
the approximation 

h 2 

ex, = £o + 777(1 - cos K z a), (3) 

where Eo = h 2 ao 2 /2m is the solution of Ec. ([2]) when 
K z ^0 and M is an effective mass [H El- Ec. (0J is the 
most commonly used expression for quasi-bidimensional 
models of superconductors (see for example (40|), but it 
is a very limited model, since it involves calculations only 
over the first energy band and with e = the ground 
state energy. In this work we use the exact solution of 
the Eq. ©. 
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A. Grand potential 

To calculate the thermodynamic properties we begin 
with the grand potential Sl(T, L 3 , which for a boson 
gas with N B particles contained in a volume V = L 3 [4l[ 
is 

fi(T, L 3 , y) = U-TS- (J,N B = ft + k BT ^ ln{l- 

cxp[-/3(e + C x {Kl + K 2 ) 1 / 2 + e K% -/*)]}, (4) 

where [/ is the internal energy, S* the entropy, \x the chem- 
ical potential, /3 = l/k B T, and the first term in the rhs 
corresponds to the K = ground state energy contribu- 
tion fi = k B T\n{l — exp[— /3(e + e — (i)]}. 

Using that ln(l — x) = — m Ec. (J4j) , and 

after some algebra we have 



f2(T, L 3 , (i) = n - k B T 



OO 

EE 

1=1 K^O 



exp{/?(/i- e )/} 



x exp{-/3(C 1 (A^ + X^) 1 ^ + 



(5) 



By substituting sums by integrals in the thermody- 
namic limit, and doing the integrals over x, y one gets 

SI (T, L 3 , fj) = k B T ln{ 1 - exp[-/3(e + e - p) } - 
I3 77 ,2 ^2 / d ^ g3{exp[-^(e A -, +eo-/x)]}- 

: i7Tj 1>1 J-00 



/3 3 (2t 



(6) 



where we have used the Bose functions g a {t) = 
Sti(^)V' CT - From ([6]) we can find the thermodynamic 
properties for a monoatomic Bose gas [34l . l35l | . 



B. Critical temperature 

To obtain the critical temperature of the boson gas 
we use the expression for the bosonic particle number, 
which we calculate summing the number of particles in 
each energy state, therefore 



Nf. 



1 



exp{/3(e +eo - ju)} - 1 



+ 



l 3 r(2) 1 

(2tt) 2 ~C 



2) 1 f°° 

YJPJ dK ^2{exp[-/3(e Kz +e - fi)}}, (7) 



where the first term of the rhs corresponds to the number 
of particles in the condensate Nq(T) and the second term 
to the ones in the excited state N e (T). 

We must notice here that from the relation for the 
Fermi energy [Il[, E F = h 2 (3ir 2 ) 2 / 3 n 2 s /3 /2m* , with E F 
the Fermi energy which corresponds to a Fermi temper- 
ature of Tp = 2290 K for the cuprate and m* = 2m e 
the effective mass of the carriers, one gets that the den- 
sity number of carriers is n s = 1.128 x 10 27 /m 3 . On the 



other hand, from the relation T /T F = 0.218 ^ one 
gets T — 499.2 K, when all the fermions in the cuprate 
are paired. Using this boson gas temperature and from 
the definition of T , one gets n B — 1.994 x 10 26 /m 3 , the 
boson density number of an ideal gas, whose value is an 
order of magnitude smaller than n s . However, as we pre- 
viously stated, analyzing the data in Fig. 2 of Ref [36| . 
and localizing the diagonal lines labeled as T = Tp and 
T = T (identified as T B ), one would expect that the 
actual quantity of superconducting carriers rib for the 
cuprate materials would be at least two orders of mag- 
nitude smaller than the n s given above. Therefore, we 
assume that only a fraction / of the maximum possible 
value n B is participating in the boson gas responsible for 
the superconductivity, so n& = fn B , and we expect this 
factor to lie in the interval f £ [0.01,0.14], according to 
what is explained in Ref. |37j . 

The BEC temperature of this boson gas in the ther- 
modynamic limit is 



r, 



1 *2 2/3 

2irn n b 



0/ mk B ((3/2)V 3 



T f 



■2/3 



(8) 



For / = 1 we recover the case where all pairable fermions 
participate in the boson gas. The corresponding thermal 
wavelenght is Ao/ = h/ ^/2-Kmk B T j = Ao// 1 / 3 . The 
quotient of the fraction of an ideal gas BEC temperature 
over the Fermi temperature of the whole is [42| 



Toy _ 2tt/ 2 / 3 

T F (6tt 2 ) 2 / 3 C(3/2) 2 / 3 



(0.218)/ 2/3 . 



(9) 



Now we use the relation for the 3D Fermi en- 
ergy Ef3D = [(37r 2 ) 2/3 /27r]£;F2D for a 3D system in 
terms of the Fermi energy for a 2D system Ef2T> — 
\m e \i 2 F2D . The constant C\ — (A/-k 2 )Ti 2 \i 2 F2D = 
[64/7r(37r 2 ) 2 / 3 ]7a 2 J B F 3 £) fc s T , where 7 = h 2 /2ma 2 k B T 
is a dimensionless constant. Introducing / in j7]), divid- 
ing by N and taking T — T c , so the chemical potential 

= £0 + eo and Nq(T c ) ~ 0, we have 



2 r(2)3xr(37r 2 ) 2 / 3 /0.436 7 T \ 1/2 



1 



/ 128 V. T F J {k B T F p c ) 2 

/ adK z g 2 {exp[-k B TF^{To/T F ){eK r - e )}}. 
Jo 



(10) 



where Ej^ t = e Az /(?i 2 /2ma 2 ), and Eq = £$/ {h 2 /2ma 2 ) 
are dimensionless energies obtained numerically from ([2]) 
for each band. 

We split the infinite integral in (ITU1) as a sum of 
integrals over the allowed energy bands, fold every 
band over the first half Brillouin zone (from to tt) 
and cut the sum at the J-th band once convergence 
has been achieved. Finally, using the numerical value 
r(2)37r(37r 2 ) 2 / 3 (0.436T /T F ) 1 /2/(i28) = 0.46532 we ar- 
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rive to 



1 = §0.46532 7 1/2 1 
V / a^ z g 2 {exp[-fc s r F /? c7 (r /r F )(e^ -e )]}, 



(11) 

which must be solved numerically. From now on we will 
be using this allowed-band-splitting method for evaluat- 
ing the infinite integrals. 




n — i i i iiiij 1 — i i i iiiij 1 — i i i iiiij 1 — i i i iiiij 1 — i i i mi 

10° 10 1 10 2 10 3 10 4 10 5 

Po 

FIG. 1: Critical temperature as a function of Po for different 
values of /. Dashed line is the experimental T c = 92.2 K for 
the cuprate. 

The experimental parameters for the cuprate 
YBa2Cu30g.92 that we use are: the critical tem- 
perature T cexp = 92.2 K [jfij : the Fermi temperature 
Tp = 2290 K @; A = 20 meV Q; the parameter 
a = 5.85 A corresponding to a = c/2 (c = 11.69 A 
the crystallographic constant), the distance between a 
copper oxide plane located at the extreme of the unitary 
cell and the Ytrium atom at the center. This is based on 
the knowledge that superconductivity occurs in the close 
vicinity parallel to the CuC>2 planes, which are two per 
unit cell. We can also calculate the thermal wavelength 
Ao = 23.575 A, and the parameters a/Ao = 0.246 and 
7 = 1.317. Finally, we obtain the magnitud of the 
jump AC/T C ~ 20 mJ/mol K 2 from the data published 
by [18j . It is well known that in cuprates there is an 
optimum oxygen dopage for which T c is higher [20j . but 
we choose the dopage for which we find more accurate 
experimental data. 

In Fig. [T] we show the critical temperature as a func- 
tion of the parameter Po for five values of /. The dashed 
line represents the experimental critical temperature for 
YBa2Cu30e.92- As we can see from this figure, there is 
only a narrow interval of values of f £ [0.001, 0.02], that 



suits the experimental condition T cexp — 92.2 K, consis- 
tent with what we expected, which in turn determines a 
set of values of Po. This means that only a small per- 
centage of the initially pairable fermions form pairs, as 
we argued above. To determine exactly both values, we 
choose a second feature of the electronic specific heat: 
the magnitude of the jump. 



C. Internal Energy 

We can derive the internal energy of the pair's gas as 
(eo + e ) , L 3 T(2) 1 



U{T,V) exp{/3(e +eo-M)}-l ' (2tt) 2 C\ /3 2 ~ 

/OO 
dK z (e + £Kjg 2 {exp[-^(e Fz + e - /i)]} 
-OO 

L 3 T(2) 2 f°° 
+2 j2^fCfW loo dK ^i e *rt-^ £K * + e o - M)]}, 

(12) 

where the first term corresponds to the particles in the 
ground state, (eo+eo)^. From the previous equation we 
subtract the ground state energy times the total number 
of pairs given by the number equation (|7|l. then multiply 
the result by / and divide it by NsksT, so we have 

(^^ £o + eo)^)^2 /27qA 
N B k B T f \T F J 



k H T, 



1 f°° 

— / adK z (e Kz - e )g2{exp[-fc B 7V/? 7 (To/T F ) ; 



(e K , + eo - fi)]} + 10.465327 1 / 2 i 

J (kbTfp) 



/OO 
adK zg3 {exv[-k B TFl3-f(To/T F )(eK z +e - /!)]}. 
-OO 

(13) 

In order to compute the internal energy and the specific 
heat it is necessary to get numerically the chemical po- 
tential and its derivative with respect to T. From the 
number equation ((TJ) and by making the considerations 
that n = fi if T < T c , and N /N b - for T > T c , fj, is 
obtained from 



1 " 7 a46532y ' 2 (OT 



adK z x 



g 2 {exp[-(k B TF(3h(T /T F )(e kz +e -£)]}, (14) 
and its derivative is given by 



T 



rf/t 
dT 



2k B T F p - ^0.46532 7 3 / 2 

J \ 1 F 



adK z x 



ln{l - exp[-fc B r F ^7(T /T F )(e ifz - fl + e )]}x 



(e Kz - n + e ) 



| .46532 7 3 / 2 



adK,x 



ln{l - exp[-k B TFpj(To/TF)(iK, - A + eo)]} 



(15) 
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III. ELECTRONIC SPECIFIC HEAT 

We consider that the electronic specific heat C e of 
the cuprate is formed by the specific heat of the gas of 
Cooper-pairs, plus the specific heat of the gas of all the 
remaining electrons that didn't form pairs. 

Like most authors do, we assume that the electronic 
specific heat is responsible for the jump, which in turn 
provides the information about the phase transition. 

In this model approach, we also consider that the spe- 
cific heat at constant volume Cy is the same as the spe- 
cific heat at constant pressure C p at least from T = 
to 200 K^is has been established for cuprates by several 
authors [2^, [44-46]. So from now on we will drop any 
subscript on this matter. 



A. Superconducting electronic specific heat 

We are able to get the superconducting electronic spe- 
cific heat C es by introducing the internal energy (TT2"j) in 
C es = [-JppU(T, L 3 )] L 3, and taking the fraction / 



C' 



1 



N b k B f \T F J ' (k B T F p) 

/oo 
adK z g 2 {exp[-k B T F ^(T /T F )(e Kz - e - ft)}} 
-oo 

x [2e Kz - e + e Q - fi + T-^} 



I / 2T N 

■27°- 46532 U ' ^ 



/OO 
adK z (e Kz - e ) 
-OO 

x ln{l - exp{-k B T F p 7 (T /T F )(e Kz + e - p)]} 

> [& , + ^ +T §, + H . 46532 ^_J^ 

/OO 
adK z g 3 {exp[~k B T F ^(T /T F )(e Kz + e - p)]}. 
-OO 

(16) 

In Figs. [2] and [3] we present the superconducting elec- 
tronic specific heat C es and C es /T as function of temper- 
ature for the gas of Cooper pairs. The height of the jump 
is reproduced by taking the values / = 0.018 (which lies 
inside the interval [0.01, 0.14] obtained in Sec. Ill B[) and 
Po = 3 x 10 5 for the dimensionless parameters of our 
model. 

The quadratic behavior of C es with temperature is 
clear from the curve presented in Fig [3j so we can ex- 
tract the slope a — 4.3 mJ/mol K 3 . It is also important 
to point out that the percentage of paired electrons that 
participate in superconductivity is among 1 and 3%, con- 
sistent with the report in 0] and in contrast to the 90 to 
95% reported in [2J|. 
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FIG. 2: Electronic specific heat as a function of T for the 
Cooper pairs gas for YBa2Cu3 06.92 using Ec. (|16[) . using 
/ = 0.018 and P = 3 x 10 5 . 
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FIG. 3: Electronic specific heat over temperature as a func- 
tion of T for the Cooper pairs gas for YBa2Cu3 06.92- 



B. Normal electronic specific heat 

We shall consider now the gas of fermions made of the 
electrons with mass m e in the spherical shell who didn't 
pair to form Cooper pairs plus the unpairable electrons, 
constituting (1 — /) of the total N electrons. The grand 
potential for an ideal fermi gas comes from (4l| 

fl(T,L 3 , Lifer) - -fc B T^ln{l + expR3(e k - M/er )]}. 

k=0 

where fif er is the chemical potential of the electron gas 
and £k = h 2 k^./2m e +h 2 ky/2m e +ek z is the energy of each 
electron free in the x — y directions and constrained by 
permeable planes in z direction. As we did in the case of 
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the boson gas, the z-component energy must come from 
the KP Eq. ^ taking the momentum of the electron 
k z and the corresponding impenetrability Pqf = Po/2. 
Replacing sums for integrals in the thermodynamic limit, 
and doing the integrals over k x , k y we have 



L 3 m e 1 



(2nf H 2 /3 2 
dkj 2 (ex.p[-(3(e kz - lifer)]), 



(18) 



where we have made use of the Fermi-Dirac functions 

f*(*) = E£i(-i) , -W 7 M- 

In order to obtain the electronic specific heat for the 
fermions C en (sometimes referred to as Cu n Q), first we 
deduce the internal energy of the fermi gas and then the 
specific heat, so doing the usual algebra we have 



= 4- 
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N fer k B ~ (1 - f) N (2irf h 2 (3 
x/ dkj 2 {exp[-/3(e kz - fi fer )]}+ 

-00 

1 L 3 m 



' (1 - /) N (2tt) 2 h 2 



dk z ln{l- 



,9flfe 



exp[-^(e/c, - ^fer)]}{2e kz - fi fer + T } 
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dk z 



(l-f)N(2nf h 2 ^ 

e kA £ k z ~ lifer + T — 



i)T 



-} 



exp[/3(e kz - li far)] + 1 



(19) 
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FIG. 4: Specific heat as a function of T for the unpaired 
electrons, using / = 0.018 and Pqf = 1.5 x 10 s . 



C. Total electronic specific heat 

The total electronic specific heat is the sum of the pre- 
vious contributions: the Cooper pairs specific heat and 
the unpaired electrons specific heat, C e — C es + C en . We 
find that the contribution of the normal state electronic 
specific heat is two orders of magnitude smaller than the 
superconducting counterpart (Cooper pairs), so in Fig. [S] 
we present the curve C e /T vs T, where the difference can 
be better appreciated, the dashed line being the normal 
electronic specific heat. 



where Nf er is the number of non-paired fermions. 

The chemical potential is calculated from the corre- 
sponding number equation as 



1 



L 3 



2 m„ 1 



(l-f)N (2tt) 2 h 2 P 
dk z ln{l + exp[-/3(e fcz - lifer)]} 



(20) 



from where Hf er and its derivative are extracted using 
numerical methods. 

In Fig. 3] we present the behavior of the normal elec- 
tronic specific heat for the unpaired electrons taking the 
values of (/, Pq) obtained in the preceding section. Notice 
that C en is linear in T until well above T c , as expected, 
and its coeficient is 7 e = 5.2 mJ/mol K 2 . 

Clearly, there is a non-zero value for 7 e , as stated in 
some reports |47| for oxygen contents x > 0.6, and it 
corresponds to the contribution of the unpaired electron, 
as suggested by Fisher, et al. in [5j. 

The analysis of the thermodynamic properties of a gas 
of fermions in a layered structure will be made elsewhere. 
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FIG. 5: Specific heat over temperature T for Cooper pairs 
plus unpaired electrons, using / = 0.018 and Po = 3 x 10 5 . 
The dashed line is the unpaired electrons contribution. 

There are some considerations we would like to do at 
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this point, based in our results. First, we do not repro- 
duce the upturn for the T < 5 K temperature region 
since the internal magnetic components and the influ- 
ence of external magnetic fields are not considered, as 
expected. Second, we have the coexistence of the two 
gases at least for temperatures under T c . Above T c paired 
fermions decouple in an unknown way, so we do not in- 
clude the decoupling mechanism in our analysis. As a 
consequence of the coexistence of the two gases, there is 
a linear temperature component in the electronic specific 
heat, which comes from (at least) the normal part of the 
system, while the quadratic temperature term is a con- 
tribution of the bosonic superconducting counterpart. In 
some articles [l2| it has been suggested that the linear 
component comes exclusively from the presence of holes, 
however, in this version of our model we can not discrim- 
inate between holes and electrons. And last, we are able 
to state that the contribution of the electronic specific 
heat to the total is around 30%, and not only 1 - 2% as 
has been repeatedly suggested [H, [23[ . 



IV. LATTICE SPECIFIC HEAT 

In this section we calculate the specific heat due to the 
lattice C;, which represents the major contribution for 
any solid. To do so, we describe and use two different 
formalisms: a simple Debye model and a phenomenolog- 
ical procedure where we take the phonon spectrum from 
inelastic neutron scattering experiments, and compare 
both results. After this, we proceed to add the result- 
ing curve to our previous electronic specific heat results. 

The total internal energy of a crystal is in general 
obtained by taking the normal vibration mode number 
G{uj)duj which lie in an interval lo, uj + cko with a fre- 
cuency u>, where G{uS) is the phonon density of states 
(PDOS). We take the energy hu> of each mode, so [ll| 



U = 



hujG{uj)duj 



{exp[huj/k B T] - 1)' 



(21) 



where we have the assumptions that the crystal is large 
enough so the sums can be substituted by integrals as 
usual. The specific heat for the lattice is 

c _ , f jhuj/k B T) 2 exp[huj/k B T}G(uj)duj 
l ~ B J (exp[Huj/k B T}-l) 2 

which is the expression from which the two methods we 
use for our analyses are derived. 



A. Debye model 

In the Debye model the solid is considered monoatomic 
and isotropic. The phonon density of states is G(uj) — 
3fc 2 /27r 2 vo = 3w 2 /27r 2 vo 3 , where vo is the sound velocity, 
which has been supposed to be the same in any direction 



(transversal or longitudinal). Introducing this PDOS in 
(f2"2"j). and after some algebra we have the Debye expres- 
sion for the lattice specific heat 



Nik B 



9s 



T 

e7> 



e °/ T x 4 exp[ X ]dx 
o (exp[x] - I) 2 



(23) 



where Q B = ^D/k B is the Debye temperature charac- 
teristic of every solid, x = fu^D /k B T \ Ni the number 
of unit cells in the solid and s = 13 is the number of 
atoms per unit cell. We must remark that in this calcu- 
lations we are assuming that the Debye temperature of 
YBa2Cu307_j; is a constant over the complete interval 
of temperatures considered. 

When one uses the Debye model for lattice calcula- 
tions, it is implicit that one is taking the harmonic oscil- 
lator approximation, so in an intent to go beyond, anhar- 
monic terms such as cubic and fourth order in T should 
be considered. However, several authors report that the 
anharmonicity does not exceed the 10 to 15% of the elec- 
tronic specific heat component in the to 300 K interval 
0, 113, so it is usually ignored. Although Debye's for- 
malism reproduces remarkably well the T 3 behavior at 
T < 5 K for cuprates, and is successful for monoatomical 
superconductors, such as Al [49j, it is not longer good 
for cuprates above 5 K, as we show in Fig. [6j where we 
present the lattice specific heat as a function temperature 
together with the experimental specific heat reported by 
Bessergenev et al for YBa2Cu306.92 [3- Even though 
the most commonly associated Debye temperature for 
this cuprate is Q B — 420 K, from now on we use the 
value ®d = 530 K in the knowledge that it provides bet- 
ter results, as will be seen later. The vertical dashed line 
in the figure indicates the critical temperature T c — 92.2 
K. 



B. Inelastic neutron scattering 

An alternative way to use series and approximations is 
to appeal to the experimental phonon density of states 
(PDOS) reported by inelastic neutron scattering (INS) 
for YBa2Cu307_a; with several doping x [Uj]. In their ex- 
periments, the authors use a non superconducting sam- 
ple as a reference, such as YBa2Cu306, which is also 
reported. 

We extract the PDOS data for YBa 2 Cu3 06.97, 
YBa2Cu307 and YBa2Cu30g from the curves given in 
[3Sj . introduce each one in Eq. (f22|) to perform the in- 
tegrals numerically observing the correct handling of the 
units, and include the curves in Fig. [6l We are assuming 
that the use of this approach already takes into account 
the anharmonic terms, at least up to the temperature 
interval considered. 

To magnify the difference among the lattice specific 
heat curves for different oxygen concentrations, we show 
in Fig. [7] the lattice specific heat over temperature for 
the same values used in Fig. [5J keeping in mind that 
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T(K) 

FIG. 6: Lattice specific heat for YBa2Cu307_ I using the 
Debye model for 9d = 530 K and inelastic neutron scattering 
(INS) for three doping x. For comparison we add the the 
total specific heat experimental data (crosses) reported for 
YBa 2 Cu 3 6 .92 by [3. 



this type of curve is the most used one for reporting elec- 
tronic specific heat. It can bee seen that the differences 
between experimental results and theoretical curves are 
more noticeable in this form. At this point we would 
like to analize three features: the difference between the 
total experimental specific heat for YBa2Cu3 06.92 and 
the lattice specific heat computed using the PDOS for 
7 — x = 6.97 at the transition point is around 30% in 
this graphic; however, the difference in the lattice specific 
heat between two adjacent superconducting dopages at 
the same point is quite small; and last, the lattice specific 
heat for 7 — x = 6 is 10% lower at T c than that for the 
superconducting counterparts. On these bases we conjec- 
ture that using the YBa 2 Cu306 compound as a reference 
for the specific heat of the lattice can be considered as a 
rough approximation. 

Above T c the lattice specific heat we obtain is still 
valid, but, as we mentioned before, our curve for the 
electronic specific heat is not strictly accurate, but only 
a guide. 

By plotting the lattice specific heat Ci/T vs T 2 , we are 
able to reproduce the T 3 term observed in some exper- 
iments [|[. As we said above, this is true for the Debye 
model, as expected, but using the INS phonon density 
of states we find a (3T 3 behavior for T < 5 K with j3 — 
0.683 mJ/mol K 4 , compared to 0.305 and 0.392 reported 
in [4| and [7fl respectively. Notice that our value of (3 
lies within the uncertainty of the INS experiment, and 
that unfortunately, no further experiments on PDOS for 
cuprates have been reported. 




T(K) 

FIG. 7: Specific heat over temperature of the phonons us- 
ing the Debye model and Inelastic Neutron Scattering for the 
same values as in Fig. [6] 

V. SPECIFIC HEAT OF YBa 2 Cu 3 06.92 

We take the electronic specific heat we calculated 
for paired and unpaired electrons with the parameters 
P = 3 x 10 5 and / = 0.018, and add them to the lattice 
specific heat from the Debye model with ©£> = 530 K 
and from the INS spectra for YBa2Cu30e.97, which we 
will be using in our calculations from now, and plot them 
in Figs. [5] and H] Together with the total specific heat 
we plot the electronic specific heat (normal plus super- 
conducting) for the purpose of emphasizing the size of its 
contribution. 




T(K) 

FIG. 8: Total specific heat for YBa2Cu3C>6.92. 

From these figures we may observe that the total height 
at the transition point for both, specific heat and spe- 
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FIG. 9: Total specific heat over temperature for 
YBa 2 Cu 3 06.92. 



cific heat over temperature, are reproduced by adding 
the three components we analized. However, we observe 
a minor difference between the experimental shape of the 
curve and ours below T c , more remarkable in Fig. [9] This 
difference becomes more notorious around 40 K, where 
we believe that the contribution of the lattice C/ needs 
a better analysis, supported by a more accurate experi- 
ment. 



VI. CONCLUSIONS 

While most authors take the experimental curves of 
the total specific heat and subtract components, we are 
able to qualitatively construct the total specific heat for 
the YBa2Cu307_a; cuprates from a simple first princi- 
ples model: the Boson-Fermion theory of superconduc- 
tivity with the layered structure of HTSC. The model 
consists in taking the Cooper pairs as a boson gas co- 
existing with a unpaired electrons (or holes) Fermi gas, 
both gases under the confinement of the layered struc- 
ture modeled by a Kronig-Penney potential in the Dirac 
comb limit in the perpendicular direction to the Cu02- 
Although no interactions among bosons and no internal 
and/or external magnetic fields are considered, the model 
reproduces qualitatively well the experimental curves of 
the total specific heat. 

A direct result is the critical temperature which de- 
pends on the anisotropy of the material, introduced 



through the planes impenetrability Pq and the separa- 
tion a between them. Since we assumed that not all 
the pairable fermions are paired, this critical tempera- 
ture also depends on the fraction of fermions that formed 
Cooper pairs. The total specific heat is calculated by in- 
cluding the specific heat coming from the bosons (super- 
conducting electronic specific heat), unpaired fermions 
(normal electronic specific heat) and the lattice. The 
resulting curve is compared to that of YBa2Cu30g.92, 
which at T c = 92.2 K has a jump AC/T C = 20 mJ/mol 
K 2 and is characterized by a linear dependence on tem- 
perature 7 e T and a quadratic one aT 2 . These last two 
features are reproducible with our model by fixing the 
parameters Pq and / with the known values of T c and 
AC/T C at T c . The values we get for the constants -f e = 
5.2 mJ/mol K 2 and a = 4.3 mJ/mol K 3 are of the order 
of the experimentally reported ones. We also show their 
correspondence with the normal electronic specific heat 
coming from the unpaired fermions, and the supercon- 
ducting term from the paired fermions, respectively. At 
the same time, these two results make plausible the as- 
sumption that not all pairable fermions in the Fermi sea 
were paired, even at temperatures near zero, and that the 
jump is a consequence of the condensation of the pairs. 
We also show that a simple Dcbye model for the lattice 
specific heat fails, as expected, but considering the results 
from INS experiments gives a better, albeit not perfect, 
approximation. A closer shape of our lattice specific heat 
to the experimental one should be obtained using data 
from a more accurate INS experiment. It can also be seen 
that the lattice specific heat shows the same temperature 
cubic behavior for T < 5 K as experiments show 

Another important result of our analysis is that the 
electronic specific heat (normal plus superconducting ) 
has a contribution of 30 — 40% of the total at the tran- 
sition temperature, and not only the 1 — 2% most au- 
thors consider, which is of the order of the specific heat 
of the unpaired fermions alone. We suggest that when 
they subtract what they consider the lattice specific heat 
from a sample, either constructed by fits or using a non- 
superonducting reference material, they might be taking 
away a representative part of the total electronic specific 
heat. 

Finally, we indirectly confirm that the upturn in the 
total specific heat at very low temperature is not a result 
of paired or unpaired fermions in the absence of a internal 
or external magnetic field. However, including magnetic 
terms should be a starting point for future study. 

We acknowledge the partial support from grants 
CONACyT 104917 and PAPIIT IN-111070 and IN- 
105011. 
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